Collective resonance modes of Josephson vortices in sandwiched stack of 
Bi 2 Sr 2 CaCu 2 08 +:1 intrinsic Josephson junctions 
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We observed splitting of the low-bias vortex-flow branch in a dense- Josephson-vortex state into 
multiple sub-branches in current-voltage characteristics of intrinsic Josephson junctions (IJJs) of 
Bi2Sr2CaCu2 08+a; single crystals in the long-junction limit. Each sub-branch corresponds to a 
plasma mode in serially coupled Josephson junctions. Splitting into low-bias linear sub-branches 
with a spread in the slopes and the inter-sub-branch mode-switching character are in good quanti- 
tative agreement with the prediction of the weak but finite inter-junction capacitive-coupling model 
incorporated with the inductive coupling. This suggests the importance of the role of the capacitive 
coupling in accurately describing the vortex dynamics in serially stacked IJJs. 

PACS numbers: 74.72.Hs, 74.50.+r, 74.78.Fk, 85.25.Cp 



Ever since the Josephson coupling was known to be es- 
tablished in Bi2Sr2CaCu208+:r (Bi-2212) single crystals 
along the adjacent superconducting CuC>2 bilayers much 
research interest has been focused on the dynamics of 
Josephson vortices in these intrinsic Josephson junctions 
(IJJs) in the long junction limit The interest 

stems mainly from that the Josephson vortex motion in 
stacked IJJs provides an ideal model to study the coupled 
nonlinear dynamic phenomena and also provides unique 
possibilities of applying stacked IJJs to high-frequency 
active devices such as sub-millimeter oscillators and mix- 
ers. 

In Bi-2212 IJJs, superconducting CuC>2 bilayers are 
much thinner than the c-axis London penetration depth, 
\ a b (~200 nm) Then, in a transverse magnetic field, 
screening supercurrents are induced in the stacked su- 
perconducting layers, which give rise to the inductive 
coupling between adjacent IJJs 5}. This coupling is 
strong enough that it leads to mutual phase locking of 
Josephson vortex motion along the c axis of stacked junc- 
tions. It has been predicted that inductively coupled and 
thus collectively moving vortices may resonate with the 
plasma oscillation modes and induce sub-branch split- 
ting of current-voltage (IV) curves near the bias voltages 
corresponding to plasma mode velocities, which was con- 
firmed by previous numerical analysis and observations 
@ . This phase-locked coherence of the plasma oscillation 
may be exploited to develop sub-millimeter oscillator de- 
vices. 

On the other hand, the thickness of CuC>2 bilayers is 
comparable to the Debye charge screening length D so 
that the local charge neutrality in the bilayers breaks 
down lj . This nonequilibrium charge variation in super- 
conducting bilayers yields adjacent IJJs are capacitively 
coupled as well. The purely capacitive coupling effect 
was treated theoretically before Recently, a model 
of the capacitive coupling incorporated with the strong 
inductive coupling has been proposed, where the capac- 
itive coupling is taken into account on an equal footing 
as the inductive one 0- In spite of the relative weak- 
ness, the capacitive coupling is predicted to lead to sig- 



nificant changes in the collective vortex dynamics. The 
model predicts an IV curve splits into linear multiple 
sub-branches with varied slopes, which corresponds to 
different plasma modes, in the low-bias region of the col- 
lective Josephson vortex flow. The spread of the sub- 
branch slopes is governed by the strength of the capaci- 
tive coupling, represented by a parameter a [=eD 2 / (st); 
e is the dielectric constant of the insulating layer, s (=0.3 
nm) and t (=1.2 nm) are the thickness of the supercon- 
ducting and the insulating layers, respectively], which 
provides a convenient means to determine the strength 
of the capacitive coupling in IJJs. 

The dynamics of Josephson vortices in stacked IJJs is 
known to be sensitive to the vortex density in the junc- 
tions 10] . The collective resonance is usually revealed 
in a high field range of H > Hd [= 4>o/2{s + t)\j) 0o is 
the flux quantum and Xj is the Josephson penetration 
depth]. In this dense- vortex state the non-Josephson- 
like emission ^lj and the Shapiro resonance steps 01 
have been observed, confirming the coherence of the 
Josephson-vortex-lattice motion over the whole stacked 
IJJs. Especially, for a higher field range of H > Hhd 
(=</>o/(s + t)Xj), the transverse inter-vortex spacing be- 
comes close to A j -"- n this range the Josephson cur- 
rent along the length of a junction distributes almost si- 
nusoidally so that the collective resonance is expected to 
be stronger. Little experimental studies, however, have 
been done for the collective motion in this highly dense- 
vortex region. 

In this study, we investigated details of multiple plasma 
modes manifested in the Josephson vortex-flow branch 
(JVFB) of stacks of IJJs in a high magnetic field above 
Hd- In all the three samples examined sub-branch split- 
ting was seen to start occurring for a field exceeding Hd- 
Especially, in the JVFB for a field above Hhd, we ob- 
tained clear hysteretic multiple sub-branches and mode 
switching between them. The linear sub-branches for dif- 
ferent collective vortex modes exhibited a spread in the 
slopes. The values of the capacitive coupling parame- 
ter a, estimated from the spread, for the three samples 
under study were within the range of the theoretical pre- 
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TABLE I: Sample parameters. N is the total number of IJJs, 
J c the critical tunneling current, (3 C [=(4J c /7r/ r ) 2 ] the Mc- 
Cumber parameter, Aj the Josephson penetration depth, and 
a the capacitive-coupling parameter. 



Sample 


stack size 


N 


Jc 


lie 


Xj 


a 


number 






(A/cm 2 ) 




W 




SSI 


17x1.5 


55 


1000 


2000 


0.31 


0.35 


SS2 


16x1.5 


25 


1500 


1300 


0.25 


0.17 


SS3 


16x1.4 


60 


1270 


2200 


0.27 


0.24 



diction of the hybrid-coupling model 9]. In addition, 
mode velocities estimated from the inter-branch switch- 
ing voltages were in reasonable quantitative agreement 
with the prediction of the model. 

Slightly overdoped Bi-2212 single crystals used in this 
study were grown by the solid-state-reaction method 0] . 
Samples were micropatterned on the crystals into the ge- 
ometry of a stack of junctions sandwiched between top 
(400 nm thick) and bottom (100 nm thick) Au elec- 
trodes [see the inset of Fig. 1(a)]. The sizes of the three 
stacks were 17x1.5 jum 2 (SSI), 16x1.5 ^m 2 (SS2), and 
16x1.4 /zm 2 (SS3). Since each superconducting bilayer 
in a stack of IJJs is much thinner than A a b, Josephson 
vortices in a usual mesa structure are supposed to be 
strongly coupled to those in the basal stack underneath 
the mesa, which may significantly distort the dynamics 
of the Josephson vortices in the mesa itself. In this study 
the basal stack was removed by using the double-side- 
cleaving technique |Tt| to overcome this problem. The 
detailed sample fabrication process is described elsewhere 
[17] . Eliminating the basal part turned out to be essential 
to obtain ideal vortex-flow characteristics of stacked IJJs. 
This sandwiched-stack geometry allowed to examine the 
collective vortex dynamics over the whole thickness of a 
stack. The magnetic field was aligned in parallel with the 
plane of junctions within the resolution of 0.01 degree to 
minimize the pinning of Josephson vortices by the for- 
mation of pancake vortices in CuC>2 bilayers [l]| • All the 
measurements were done at 4.2 K in a two-terminal con- 
figuration [Fig. 1(a)]. The combined contact resistance 
for bottom and top interfaces (~2 kf2 for SSI and SS2, 
and ~600 f2 for SS3) was subtracted numerically. 

The sample parameters are listed in Table I. The num- 
ber of IJJs in a sandwiched stack, N, was estimated 
from the number of zero-field quasiparticle curves (not 
shown) . The McCumber parameter (3 C was obtained from 
the zero-field return current I r [as denoted by arrows in 
Fig. 1(a) for decreasing bias]. Fig. 1 shows the JVFB of 
sample SSI with a single up-down bias sweep. The values 
of Hd and Hhd for SSI estimated with Aj were 2.2 and 4.4 
T, respectively. Thus, the data in Figs. 1(a) and 1(b) cor- 
respond approximately to the two dense-vortex ranges, 
H ~ Hd and H ~ Hhd, respectively. The average slope, 
AV/AI, of the JVFB keeps increasing in proportion with 
the magnetic field intensity. A few voltage jumps are seen 
in the JVFB for fields beyond ~2.6 T in the bias range 
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FIG. 1: (color). Josephson vortex-flow branches of SSI in 
field ranges of (a) H ~ Hd and (b) H ~ Hhd (the 4-T curve is 
shifted upward by 1 jj, A). Neighboring sub-branches are col- 
ored alternately for clarity. The arrows in (a) indicate the cur- 
rent sweep direction. The arrows in (b) indicate the bound- 
aries between the vortex-flow and the quasiparicle regions, 
defining the parameter . Inset of (a): the measurement 
configuration. Inset of (b): low-bias vortex-flow region at 5.2 
T 



above ~40 mV. This unusual feature, taking place in the 
vortex-flow state for a bias below the quasiparticle-state 
return current l r in each field, is caused by switching of 
coherent vortex-lattice motion between different plasma 
modes. Coherently moving vortex lattice remains in a 
certain transverse plasma mode as long as the vortex- 
lattice velocity is smaller than the propagation mode ve- 
locity of the plasma oscillation. Exceeding the mode ve- 
locity for a higher bias current, however, the resonating 
dynamic state of the Josephson vortex lattice suddenly 
switches to the adjacent plasma mode with a higher prop- 
agation velocity, giving rise to a voltage jump observed. 
The onset field of the observed voltage jumps almost co- 
incides with the value of Hd- 

Clearer voltage jumps in the JVFB are seen for a higher 
vortex density in the field range of H > Hhd as in Fig. 
1(b), as the resonant coupling between the vortex-lattice 
motion and the plasma oscillation is strengthened. In this 
field range the hysteresis in the quasi-particle branches 
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for V > V£f (denoted by the arrows) in Fig. 1(b) is 
almost completely suppressed while the voltage jumps 
become more evident in JVFB. The inset of Fig. 1(b) 
shows the distinct voltage jumps in the lower bias region 
for 5.2 T. One notes that, in Fig. 1(b), all the voltage 
jumps take place in the low-bias flux-flow region bounded 
by V£f . The voltage interval between the neighboring 
jumps tends to increase for higher biases, which is con- 
sistent with the prediction of both the purely inductive 
and the inductive-capacitive hybrid coupling models. 

Fig. 2(a) displays the JVFB of the sample SS2 (with 
7V=25) taken for =5.8 T. For the measurement the bias 
current was repeatedly swept up and down so that de- 
tailed hysteretic feature of the sub-branches were traced 
out. These multiple-sweep data clearly indicate that the 
voltage jumps observed in the single-sweep data as seen 
in Fig. 1(b) were indeed caused by the inter-sub- branch 
mode switching. Since the Hhd for this sample is 5.4 T 
the data in Fig. 2 correspond to the highly dense vortex 
state of H > Hhd, expectedly with many resonant vor- 
tex states as in the sample SSI in the similar field range. 
About 22 multiple sub-branches were obtained in SS2, 
which was close to the number of the IJJs iV=25 in the 
stack. Details of the vortex- flow region of 5.2 T curve of 
SSI, as shown in Fig. 1(b), are also illustrated in Fig. 
2(b) for a quantitative analysis. 

In Figs. 2(a) and 2(b), all the discernible sub-branches 
exhibit linear IV characteristics in the low-bias vortex- 
flow region for V < V^ , the occurrence of which is 
in agreement with the prediction of the induct ive-and- 
capacitive hybrid coupling model. By contrast, the 
purely inductive coupling yields sub-branch splitting near 
the voltages corresponding to the plasma mode veloci- 
ties. As denoted by the grey region in Fig. 2 the linear 
sub-branches, when extrapolated to the low-bias region, 
converges to a single intercept point with a positive value 
Ip on the current axis. It confirms the vortex-flow nature 
of the sub-branches with finite vortex pinning, that may 
have been caused by any defects in the IJJs or by any 
pancake vortices formed in Cu02 layers by a misaligned 
field component. The pinning current I p was 2.07 fiA and 
0.64 fiA in SSI and SS2, respectively. For an analysis of 
the finite spread in the slopes of the linear sub-branches 
we adopt fitting to the hybrid coupling model 0, while 
incorporating the effect of the finite vortex pinning (rep- 
resented by I p ) and the change of the vortex density with 
applied magnetic fields (H ext ) as 



V = N(I-I P )[^ + ^^] 

A n l c 



(1) 



Here, A n = l + 2a[l — cos(nir / (N + 1))], I c is the critical 
Josephson current at a given applied field, and k is a 
constant. Since samples are in the quasiparticle state for 
V > V^ we define the current corresponding to V£* in a 
given field as the critical current I c in Eq. (1). Also, the 
vortex velocity for / > I p is assumed to increase linearly 
with the bias current. 

Without the capacitive-coupling effect (a=0) the low- 
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FIG. 2: (color). Josephson vortex-flow sub-branches (a) for 
SS2 at 5.8 T and (b) for SSI at 5.2 T with the spread of linear 
sub-branches denoted by the gray region. Insets: details of 
the sub-branches in the low-bias range. 



bias sub-branches with spreaded slopes, as exhibited by 
the grey regions in Figs. 2(a) and 2(b) for SSI and 
SS2, respectively, are supposed to reduce to single lines 
(pointed by arrows) , where the value of k determines the 
slope of the lines. Increasing a the spread of sub-branch 
slopes increases from the a=0 lines to the left. In the 
very low bias region the JVFB in both samples becomes 
chaotic. 

The parameter a was determined from the best fit of 
the spread to Eq. (1), where the normal-state resistance 
R n was fixed separately following the fitting procedure in 
Ref. [H. The best-fit values of a for SSI at 5.2 T, SS2 
at 5.8 T, and SS3 at 4.15 T were 0.35, 0.17, and 0.24, 
respectively. These values are in the range of theoretical 
estimate of the hybrid coupling model 0, 0.1 ~ 0.4, 
and in reasonable agreement with the recently observed 
values of 0.36 - 0.44 in SmLai^S^CuO^a 20]. The 
corresponding values of the k for SSI, SS2, and SS3 were 
0.32, 0.14, and 0.38 mV/T, respectively. 

Fig. 3 displays the mode-switching voltages (filled cir- 
cles) and boundary voltage V£* (filled square) for SS2 
as a function of the mode index n. Out of 25 expected 
modes corresponding to N the three lowest-index modes 
were not observed, however, presumably because they 
were located beyond V£* . The solid line is the trace 
calculated with a=0.17 and the best fit value for the 
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FIG. 3: Mode-switching voltages (filled circles) and 
(filled square) as a function of mode index n for SS2. The 
solid line is the best fit to the hybrid coupling model. The 
dotted line is obtained for the inductive coupling only, with 
the same parameters as used for the solid line. Inset: the best 
fit to the purely inductive coupling model. 



Swihart velocity cq =1.05x10 5 m/s. The resonance volt- 
age of the nth mode is obtained from the relation 
V n = Nc n H ext (t + s). The mode velocity of the trans- 
verse plasma oscillation (c„) is expected to be enhanced 
by a factor of \J A n from that (c°) of the purely inductive 
coupling mdoel (A n = 1) by the nonvanishing capacitive- 
coupling effect |j| as in the relation 

Cn = C \ - . A " — rr = C° n y/A^ , (2) 

y 1 — cos(7m/ (TV + 1)) 

where the mode index n runs from 1 to N . The dotted 
line corresponds to the case of purely inductive coupling 
(cv=0) with the same Swihart velocity as used for the 
solid line. One sees that the enhancement of the mode ve- 
locities due to nonvanishing capacitive coupling becomes 
more evident for high mode indices, giving a maximum 
enhancement for cn as cn — [(l + 4a)] 1 / 2 . The inset of 



Fig. 3 illustrates the best fit of the mode voltages to the 
purely inductive-coupling model, with the best-fit value 
of Co to be 1.35xl0 5 m/s. In this case, mode voltages 
for low-index modes are seen to fall significantly below 
the model expectation. This trend was also observed 
previously in the motion of microwave-induced Joseph- 
son vortices, although the vortex density was much lower 
than in this study [l5|. 

The Swihart velocity for SS2, estimated using the re- 
sistively shunted junction model, was cq = 2-7r/Aj=8.7x 
10 4 m/s. In the calculation we used the Josephson 
plasma frequency f p i (—\/ I c /2tt^> Cj) of 55 GHz, es- 
timated from the critical current / c =0.358 mA and the 
junction capacitance Cj=9 pF determined from (3 C and 
R n . The value of cq is close to the one obtained for the 
hybrid-coupling model but tends to deviate more from 
the predication of the inductive-coupling model. 

In summary, the observed characteristics of the sub- 
branch splitting in the vortex-flow region were in reason- 
able quantitative agreement with the inductive coupling 
incorporated with the capacitive coupling. The multi- 
ple sub-branches in JVFB were shown to be caused by 
the inter-junction capacitive coupling and the value of 
the coupling parameter a estimated from the spread of 
multiple branches was in good agreement with the the- 
oretical prediction of the hybrid coupling as well as the 
previous observation in optical measurement [20|. The 
quantitative identification of the multiple sub-branches 
in this study confirmed the significance of the capacitive 
coupling in addition to the inductive coupling in accu- 
rately describing the dynamics of Josephson vortices in 
serially stacked IJJs. For further detailed analysis of the 
Josephson vortex dynamics consideration of the charge- 
hole imbalance coupling [2l| and the dissipation effect of 
quasiparticles in the ab plane [22| may be required. 
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